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The possibility of a quartet superfluidity in fermionic systems is studied as a new aspect 
of atomic gas at ultralow temperatures. The fourfold degeneracy of a hyperfine state and 
moderate coupling is indispensable for the quartet superfluidity to occur. Possible supercon- 
ductivity with quartet condensation in electron systems is discussed. 
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The superfluidity is so far known to be sustained by the Bose-Einstein condensate or the 
Cooper-pair condensate. The former is realized in liquid 4 He and some atomic gas of an alkali 
element, 1-4 while the latter is realized in liquid 3 He and a variety of superconductors. Recently, 
it has been found that Cooper pairs are formed in the fermionic atomic gas of alkali element 5 
and the crossover to the Bose-Einstein condensation of a diatomic molecule is possible with 
the help of Feshbach resonance. 6 

In principle, there exists the other possibility that superfluidity is sustained by a conden- 
sate based on four fermions (quartet) as in the a-particle correlation in a light nucleus. 7 An 
a-particle consists of two protons and two neutrons which have approximately 4-fold degener- 
acy corresponding to 2 x 2 degeneracy of the real spin and isotopic spin states. Such 4-fold de- 
generacy is possible in a fermionic atomic gas with total angular momentum F = S + I = 3/2, 
S and I being the electron spin and nuclear spin, respectively. For example, a Be atom has 
F = 3/2 with S = and / = 3/2. 

The purpose of this Letter is to assess the possibility of quartet superfuidity by analyzing 
an extended Cooper problem, 8 i.e., a four-fermion problem outside a rigid Fermi surface. This 
is regarded as a natural extension of the Cooper problem to the present case. We expect that 
the same tendency as in the "Cooper problem" is taken over to the many-body problem in 
which the deformation of the Fermi surface should be taken into account self-consistently. 
One might think that there is another way of extending the Cooper problem to the present 
case, i.e., by analyzing a stability condition of the BCS state against an addition of four extra 
fermions to the BCS state. However, we do not take such an approach here because not only 
mathematical treatments but also the concept of a theoretical framework seems to be much 
more complicated than ours. The result of this study might also give us a clue to the possibility 
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of quartet superfluidity in the electronic system with a nearly degenerate Fermi surface. 

A theorem by Nagaoka and Usui 9 states that the ground-state wave function of a many- 
particle system is symmetric with respect to an interchange of position between any pair of 
particles if one could neglect the possible spin dependence. Therefore, the ground state of four 
fermions with "spin" F = 3/2 should be 

* Q (1, 2, 3, 4) =$(r!,r 2 ,r 3 , r 4 )x(l, 2, 3, 4), (1) 

where $ is a fully symmetric function with respect to the interchange of f. 
wave function for the spin is given by the fully antisymmetric one: 

a(l) f3{l) 7 (1) 5(1) 

\_ a(2) (3(2) 7 (2) 5(2) 

Vi! a(3) (3(3) 7 (3) 5(3) 

a(4) (3(4) 7 (4) 5(4) 

Here, a, (3, 7, and 5 denote the "spin" states F z = 3/2, 1/2, —1/2, and —3/2, respectively. 
Due to the Pauli principle, the orbital part 3>(ri, ra, 7*3, rk) is a fully symmetric function of rj. 
State (1) represents that of a tetra-atomic molecule (quartet) if the bound-state solution is 
possible, and is more stable than that of two diatomic molecules whose orbital wave function 
is not fully symmetric. 

Now, we discuss a four-body problem outside the rigid Fermi surface, namely, the Cooper 
problem. Here, we compare the energy of the quartet and two Cooper pairs. Because it is 
difficult in general to obtain the exact solution of a four-body problem, we try to solve the 
problem variationally. The variational solution for <3? is assumed to be in the Hartree form as 



X(l,2,3,4) 



(2) 



$(ri,f 2 ,f 3 ,f4) = Y[f(n-R), 



(3) 



i=l 



where R = Yli=i is the center of gravity. The single-particle state f(r) is expressed as 

m= E 4- eife> ( 4 ) 
ifci>fe F 

It is noted here that only the states outside the Fermi surface \k\ > kp are available. 
The energy for state (3) is given as 



Eq = [■■■ fdr 1 ---dnl[f(7 

J J i=i 



R) 



i=l 



i>j 



i=i 



= 4 J dff(r) [-^ + || df>V(r- 7)\f(r 



f(r), 



(5) 



where V(fi — fj) denotes the two-body interaction between fermions. Optimization of Eq with 
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respect to f(r) leads to the equation 

h 2 v 2 



2m 



+ 



3 J dfV(f- 01/(0 



/(r) = A/(r), 



(6) 



where A is the eigenvalue. This equation should be solved under the normalization constraint 



/ 



dr[/(r)] 5 



1. 



The resultant expression of optimized energy is given by 

h 2 v 2 



A 



J dr/(r) f- 



2m 



+ 4 



y df/(r) 



h 2 v 2 

2m 



(7) 



(8) 



In order to discuss the "Cooper-type problem", let us introduce the ^-representation for 
the interactions as 



and 



U(r) = 3 J drV(f- Ol/(0! 2 = E^ 9 "' 



where [7^ is expressed in terms of V$ and as 

Uq = 3V^^ ] fq-pfp- 

p 

Then, the effective Schrodinger equation for the 'Cooper problem" is given as 



(9) 



(10) 



(11) 



(12) 



where e = h 2 k 2 /2m. It is noted that the wave vector in the wave function is restricted in 
a way such that /cf < |fe|, A;f being the Fermi wave number. 

In order to solve eq. (12), following BCS, we introduce a model attractive interaction as 
follows: 



\-V, (e F < e fc ,e fc / < e c ); 
[ 0, (otherwise). 

Then, the second term of eq. (12) is expressed as 



E u k-«fk> = ^l k kc dk'k' 2 f(k')J^dtu(\k-k'\ 

= dk'U*(k,k')k' 2 f(k f ), 

Jk F 



(13) 



(14) 



where \k-k'\ = (k 2 + k' 2 - 2kk't) 1 / 2 , i.e., t = (k ■ k')/kk', and the kernel U* is defined as 

1 



U*(k,k') 



4vr 2 



dtU(\k-k'\). 



(15) 
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Here, the function U(\k — k'\), (11), is expressed as 

U(\k- k'\) = SV U _~ f_&t' j\pp 2 f(p)f{\k- k'-p\), (16) 

where t' = (p- (k — k'))/p\k — k'\. 

Finally, by introducing the function F(k) = kf(k), eigenvalue equation (12) is transformed 
to the following form: 

e k F{k) + j ° dk'kU*{k, k')k'F(k') = XF(k). (17) 

Now, the kernel kU*(k,k')k' is symmetric with respect to the interchange k <-> k'. Then, the 
eigenvalue problem of eq. (17) is easily solved numerically if the "interaction kernel" U*(k, k') 
is given. We solve numerically the self-consistent set of equations, (15)~(17). Then, the energy 
of the quartet (eq. (8)) is given as 

(E Q ) miQ = 2^X + j" |* efc |F(A:)| 2 ) . (18) 

This energy should be compared with that in the original Cooper problem. The energy of 
one Cooper pair Eq is given by 



— I 



fee h.2 

dk -. (19) 

The phase diagram in the pYV-{k c /k-p) plane, pf = mk-p /2TT 2 t^ , is given in Fig. 1 in which 
equi-energy surface of (-Eq — 2Eo)/e-p is shown. It is seen that the quartet state is stabilized 
in the intermediate and strong coupling region with k c > 3/cf- This is consistent with the 
Nagaoka-Usui theorem which indicates the quartet state to be more stable than the two- 
Cooper-pair state in the dilute limit, i.e., without the Fermi see. It is noted that there exists 
a threshold for eq. (17) to have a bound-state solution, while there exists no threshold for 
the original Cooper problem (eq. (19)). Therefore, in the limit of weak coupling, the fully 
symmetric solution in the orbital space will not be the ground state of the "Cooper problem" , 
in contrast to the 4-body problem in vacuum. It is noted that the quantitativeness of the phase 
boundary in Fig. 1 may not be very good because we adopted the Hartree approximation 
to solve the 4-body problem and neglected the interaction between two Cooper pairs as in 
the BCS approximation. Nevertheless, this gives us a guideline to understand the crossover 
between the Cooper-pair and the quartet condensations. 

Although we adopted here the model interaction (13), it is not so difficult to discuss on 
the basis of a more realistic interaction similar to the 12-6 Lennard-Jones potential, as we 
will discuss elsewhere. Since the diatomic bound state is barely possible for the 4 He system, 10 
the weakest interaction and lightest mass, it is possible for a neutral atom with a stronger 
attraction and a heavier mass, such as Be atom, to form the diatomic bound state. Then, 
due to the Nagaoka-Usui theorem, the quartet bound state is also expected to be possible. In 
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addition, the tuning of attractive interaction may be possible with the use of a mechanism 
similar to the Feshbach resonance for two-body interaction. 11-13 This would make it possible 
to search a crossover phenomenon between the Bose-Einstein condensation of tetra-atomic 
molecules and the quartet Fermi superfluidity. 

In order to discuss the relative stability between the Cooper-pair and quartet condensed 
state, we of course must develop a self-consistent treatment that relaxes the rigidity of the 
Fermi surface. Let us briefly sketch how to describe the many-body quartet condensed state. 
A creation operator of a 4-atom molecule with zero total momentum is given as 

b] = E • • • E +^2 + k 3 + kMh,k 2 , h, ^)4 ia 4 2/3 4 37 4 45 . (20) 

ki fc4 

It is natural to assume that the wave function is given by a product of a one-particle state 
as in eq. (3): 

4 

Hh,k 2 ,k 3 ,k 4 ) = H<p(k l ), (21) 
i=i 
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Then the many-body quartet condensed state |^qc) may be written in a manner analogous 
to the BCS state as 

I*QC) = Vn/4 exp (^Jjb^J |vac), (22) 

where N is the total number of fermions, the projection operator 7-V/4 works to project out 
the iV/4-quartet state, and |vac) denotes the vacuum. State (22) is equivalent to 

{fei~fc 4 } i=i ' 

where the factor S(ki + ki + £3 + £4) has been abbreviated in the product Y\t=n an d 
4>{k) = v(k)/u(k) and the normalization condition |u(/e)| 2 + |v(A;)| 2 = 1 should be fulfilled. The 
variational functions u(k) and v(k) should be determined so as to minimize the free energy, 
as in the BCS theory. 14 

Finally, let us discuss the possibility of observing quartet superconductivity in solids. The 
most promising case is in metals with two almost degenerate bands, e.g., the a- and 6-bands, 
and an extremely low carrier density. In such a system, quasi-4-fold degeneracy is maintained 
for each k-point by spin a = ± and band index m = a,b. Diluteness of carriers makes it 
possible for the plasmon mechanism to work, 15,16 and the energy range of attraction huj p i 
extends over the Fermi energy ep: 



» 1, (24) 



where ro and ae denote the mean distance between electrons and the Bohr radius, respectively. 

In conclusion, the possibility of quartet Fermi superfluidity and superconductivity have 
been investigated on a model interaction of BCS type. 
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Note added: 

After the present paper had been accepted for publication, we became aware of the pre- 
vious works which discussed the physics similar to the present one. Transition temperature 
of the quartet condesation in nuclear matter has been discussed in G. Rdpke, A. Schnell, P. 
Schuck and P. Nozieres: Phys. Rev. Lett. 80 (1998) 3177. Model system in the large spin limit 
has been analyzed in A. S. Stepanenko and J. M. F. Gunn: cond-mat/9901317. The problem of 
one-dimensional models has also been discussed in P. Schlottmann: J. Phys.: Condens. Matter 
6 (1994) 1359; and C. J. Wu: cond-mat/0409247. We thank Congjun Wu who directed our 
attention to those works. 
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